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(3 1 INTRODUCTION 

; The supernova remnant G21. 5-0.9 (SNR 021.5-00 . 9) ha s 
, been known for over 30 years (e.g., lAhenhoff et aLlll970t ). 
5^ and has long bee n classified as a filled-c entre or a Crab- 
\ like remnant (e.g.. Iwilson fc Weileij[l976l ). Such supernova 
remnants are powered by a central pulsar rather than by 
the interaction of the expanding ejecta shell with its sur- 
roundings. We will refer to them using the term "pulsar 
wind nebula" or PWN. G21.5— 0.9 is a bright, centrally con- 
densed radio and X-ray source, with a diameter of ~l'. A 
larger, low sur face brightness h alo, of diam eter ~2.5', is seen 
in the X-rav |Bocchino et al.|[20 05: Mat heson fc Safi-Harbl 
I2OO5I : ISafi-Harb et al.ll200ll: ISlane et al..i200ol). The distance 



ABSTRACT 

We report on new 5-GHz VLA radio observations of the pulsar-powered supernova 
remnant G21.5— 0.9. These observations have allowed us to make a high-quality radio 
image of this remnant with a resolution of ~0.7". It has a filamentary structure similar 
to that seen in the Crab Nebula. Radio structure suggestive of the torus seen around 
the Crab pulsar is tentatively identified. We also compared the new image with one 
taken ~15 yr earlier at 1.5 GHz, both to find the expansion speed of the remnant 
and to make a spectral index image. Between 1991 and 2006, we find that the average 
expansion rate of the remnant is 0.11 ± 0.02% yr~^, corresponding, for a distance 
of 5 kpc, to a speed of 910 ± 160 km s^^ wrt. the centre of the nebula. Assuming 
undecelerated expansion, this expansion speed implies that the age of G21.5— 0.9 is 



yr, which makes PSR J1833— 1034 one of the youngest, if not the youngest. 



870 

known pulsars in the Galaxy. 
Key words: ISM: supernova remnants 



to G2 1. 5-0.9 is ~5 kpc l|Camilo et al.ll2006l : ]Davelaar et al.l 
1 19861 ). Its radio and X-ray luminosities are ~10% and ~1% 
of those, respectively, of the Crab Nebula. 

However, despite numerous searches, no pulsar was 
seen in G21.5— 0.9 until very recently, when t wo teams an- 
nounced the discovery o f PSR J1833-1034 (|Gupta et al.l 
I2OO5I : ICamilo et al.|[2006l ). The pulsar is very faint, with a 
pulse-averaged fiux density of ~ 70 /iJy at 1.4 GHz. It has 
a period of 61.8 ms and a period derivative of P = 2.02 x 
10~^^, giving it a characteristic age of 4800 yr. Although 
no other relatively direct measurements of G21.5— 0.9's age 
have been made, various ar guments have led other authors to 
suggest an age of ~1000 yr ( Bock et al.ll200ll : lBocchino et al.l 



l2005l ; ICamilo" et al.|[200^ . Despite its faint pulsed emission, 
the pulsar is in fact very energetic, having a spindown lu- 
minosity oi E = 3.3 X 10^^ erg s~^, which, in our Galaxy, is 
second only to that of the Crab pulsar. 

The powerful winds of young pulsars with high spin- 
down luminosities have been shown to excite complex 
structures, visible fro m the X-ray to the radio (e.g., 
ICaensler fc SlarL3l2006l ). In particular, in the Crab Nebula, 
narrow emission features, called "wisps" are seen. The wisps 
are thought to be associated with the shock in the pulsar 
wind, and with a torus seen in the X-ray. They move rapidly 
outward and are seen in the radio, optical and X-ray (e.g., 
iBietenholz et al.ll2001al . l2004 and references therein). 3C 58, 
which is also the wind nebula of a young, energetic pulsar, 
exhibits a narrow radio feature which is probably associated 
with the pul sar outfiow, although it is not as mobi le as the 
Crab wisps (|Frail fc Moffettlll993l : lBietenholdl2006l ). 

While numerous radio images of G21.5— 0.9, at 
frequencies from 330 MHz to 94 GHz, have been 



published (|Becker fc Kundu 


Il976l: IWilson fc Weiled 


19761: Becker fc Szvmkowiak 


198ll: Morsi fc Reich Il987: 


Fiirst et al.l 1 19881: iKassimI 1 19921: iBock et al.l I2OOID. none 



of them had a resolution higher than 4", which is not 
sufficient to see such wisps and other radio features 



detail. iFrail fc MoffettI (|l993l ) observed G21.5-0.9 with the 
VLA specifically to search for radio emission associated 
with the (then still unknown) pulsar, but found none. Now 
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the pulsar location is known, it seemed well worthwhile to 
revisit G21.5— 0.9 with deeper and higher-resolution radio 
observations. We therefore observed G2 1.5— 0.9 using the 
NRAC0 Very Large Array (VLA) with the goals of making 
a sub-arcsecond resolution radio image, searching for fea- 
tures near the pulsar, and, by comparing our observations 
to older ones, determining the expansion rate, age, and 
spectral index distribution of the synchrotron nebula. We 
describe the observations in § [51 present the new image in 
§ |3l determine the expansion rate and form a spectral index 
image in § |4l and discuss our findings in § [5] 



2 OBSERVATIONS AND DATA REDUCTION 

We observed G21.5— 0.9 at 5 GHz, using the A and B array 
configurations of the VLA. All our observations were phase- 
referenced to the radio source PMN J1832— 1035. The details 
of the observations are given in Table[TJ which also describes 
the two archival data sets that we use in this paper. The flux 
density calibration for all observing sessions was done by 
observing 3C 286. The data reduction was carried out using 
standard procedures from NRAO's AIPS software package, 
which were also used to re-reduce the archival data. 

To make the best possible image, we used maximum 
entropy deconvolution, which is better suited to extended 
objects like G21.5— 0.9 than the more common CLEAN de- 
convolution. However, since the diameter of the PWN is ^1', 
even the shortest spacings with the B-array recover only 
around half the total flux density. Therefore the large-scale 
structure may not be well recovered without some additional 
short-spacing information. We used an archival 5 GHz VLA 
data set taken in the D configuration on 30 Nov 1985 to 
supply additional short spacing information. 



in rad io images o f the Crab nebula (iBietenholz et al. Il2001al . 
|2004) and 3C 58 (|Bietenholzll2"00^ : IBietenholz et al.ll20oTbh . 
is clearly visible. 

The remnant has a double-lobed structure, with the axis 
of symmetry running approximately from the northwest to 
the southeas t. This structure has also b e en noted in the 
X-ray (e.g., iMatheson fc Safi-Harbl l2005l : ISafi-Harb et~al] 
l200li ). Unlike in the X-ray, however, there is no bright cen- 
tral condensation in the radio, and the northwest lobe is 
somewhat brighter than the southeast one. 

The position of the pul sar is known to within < 0.5" 
from Chandra X-ray images (jCamilo et al.| [2006fl. Although 
the pulsar is too faint to discern in our radio images, we 
can accurately locate it, since our observations were phase- 
referenced to the radio calibrator PMN J1832-1035. We 
estimate the uncertainty of our astrometric calibration is 
also about 0.5". The pulsar position is indicated in Fig. [1] 
The pulsar is displaced by ~4" southeast of the centroid 
of the remnant. Loop-like filaments are seen, with the most 
prominent being to the northwest. 

We show a detail of the region near the pulsar, made 
with a slightly higher resolution (AIPS robustness parame- 
ter = -1) of 0.55" X 0.37" at p. a. 7° in Fig.[l The brightest 
region is ~10" north-northwest of the pulsar. The pulsar 
itself is in a region of somewhat lower surface brightness 
of diameter ~8", which is open to the southeast. There is 
possibly an elliptical structure, incomplete to the southeast, 
oriented at p. a. 45° and of major axis diameter ~8" (cor- 
responding to 6 X 10^^ cm at 5 kpc), although a conclusive 
identification is not possible at this combination of signal-to- 
noise and resolution. If the structure is real, then its minor 
axis is oriented along the axis of symmetry of the nebula 
(i.e., northwest to southeast). 



3 HIGH RESOLUTION IMAGE 

In Figure[T]we show the full radio-image of G21.5, obtained 
by combining the A, B, and D array observations at 5GHa 
and using a AI PS' implementation of the Briggs robust- 
ness parameter (|Briggs et al.l Il999l ) of 1.5, which achieves 
a resolution somewhat better than available through nat- 
ural weighting with only slightly increased noise. The full- 
width at half-maximum (FWHM) of the elliptical Gaussian 
restoring beam was 0.82" x 0.53" at p.a. 10°. The total 
flux density recovered was 6.7 ± 0.3 Jy, consistent with the 
earlier total flux density measurements at cm wavelengths 
l| Wilson fc Weileij|l976l : ISalter et al.lll989l ). The off-source 
rms background brightness was 19 /iJy bm"^ This image 
has notably higher resolution and lower rms background 
than any previously published radio images of G21.5— 0.9. 
Prominent filamentary structure, reminiscent of that seen 

^ The National Radio Astronomy Observatory, NRAO, is a facil- 
ity of the National Science Foundation operated under coopera- 
tive agreement by Associated Universities, Inc. 
^ In anticipation of our result on the expansion, described in ii l4.1l 
below, the D-configuration data, taken 21 years before the A and 
B-configuration data, was scaled up by 2.3%, based on our mea- 
surement of the expansion speed. The effect on the image of this 
small correction to the D-configuration data is minimal. 



4 EXPANSION SPEED AND SPECTRAL 
INDEX 

We set out to both determine G21.5— 0.9's rate of expansion 
and to map its radio spectral index by comparing our 5 GHz 
image from 2006 to a 1.5 GHz one ma de from data taken in 
1991 in the A array configuration (see lFraiFfc Moffett 199^, 
for the original results from the 1991 data). Although an 
archival image at our observing frequency would have been 
preferable for determining the expansion, no older 5-GHz 
image of sufficient sensitivity and resolution was available. 
Since the two images differ in both epoch and frequency, 
the expansion and spectral index cannot be determined in- 
dependently of each other. However, for a source such with 
well-resolved, small-scale structure such as G21.5— 0.9, and 
for a range in frequency of only three, it is highly unlikely 
that variation in the spectral index would cause the source 
to appear larger or smaller at different frequencies. Further- 
more, earl ier results on G21. 5-0.9 l|Morsi fc Reicbl Il987l : 
iFiirst et"al . 1988) suggest that G21.5— 0.9's radio spectral 
index is relatively uniform across the nebula. Large varia- 
tions in spectral index are not in fact expected in PWNe, if 

We note that the position of the pulsar itself has not been ac- 
curately measured, but that the position of the bright, although 
un-pulsed, cent ral condensation in X -rays is identified as the pul- 
sar position bv lCamilo et al! I ||2006|') . 
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Table 1. Observing Sessions for G21.5— 0.9 



Date 


Frequeney (GHz) 


Array Configuration 


Observing Time'' (hrs) 


2006 Mar 7 


4.93, 4.56 


A 


6 


2006 Sep 9 


4.93, 4.56 


B 


2 


1991 Jul 3<=<* 


1.46, 1.51 


B 


6 


1985 Nov 30<= 


4.89, 4.84 


D 


1.5 



^ The sky frequencies of the two intermediate frequencies (IF) used. In each case the 
bandwidth per IF was 50 MHz. 

^ The approximate total observing time, including calibrator observations; note that for 
1985 Nov 30, only "^10 min were spent on-source for G21.5--0.9. 

Data from the NRAO data archive 

Original results published in lFrail fc MoflettI ^^). 




RIGHT ASCENSION (J2000} 

Figure 1. A radio image of G21.5— 0.9 at 4.75 GHz taken on 2006 March 3. The FWHM of the elliptical Gaussian restoring beam was 
0.82 X 0.53" at p. a. 10°. Contours are drawn at —4, 4, 8 and 16% of the peak brightness of 1.69 mjy bm~' and the rms background 
brightness was 19 f iiy bm"'^. The grey scale is labelled in mJy bm"'^. The white cross marks the position of the pulsar, taken from 
ICamilo et al.l ||2006| ). and known to an accuracy of < 0.5". The image was made from A and B array configuration data taken in 2006, 
with the addition of a small amount of archival D array configuration data from 1985 (see text, § [2]l. In the electronic edition, the 
accompanying animation shows alternating images from 2006 and 1991, visually showing the expansion over the 15-year period. The two 
images are on the same scale, and the larger and smoother of the two is that from 2006, with the smoother appearance being due to the 
higher signal-to- noise (animated gif file available from jhttp : //lanl . arxlv . org | by downloading "source" from "other formats"). 
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Figure 2. A detail image of tlie region near the pulsar, made 
from the same data but at a slightly higher resolution than that 
of Figure [T] The axes show the RA and decl. offsets from the 
pulsar position, which is marked with a cross. The FWHM of the 
elliptical Gaussian restoring beam was 0.55 X 0.37" at p. a. 7°. 
The grey scale is labelled in /^Jy bm~^. 

the remarkably uniformity of spectral index over each of the 
two most pr ominent remnants, the Crab N ebula and 3C 58, 
is any guide l|Bietenholz et al" ]|l997l . l2001bl ). Therefore, since 
the expected variations in spectral index are small, and very 
unlikely to mimic a frequency-dependent change size, we 
proceed to compare the 2001 1.5-GHz and 2006 5-GHz im- 
ages to determine the expansion rate. Once the expansion 
rate is known, we can "expand" the 1991 1.5-GHz image in 
order to make a spectral-index image. 

4.1 Expansion 

We used th e same appr o ach to det ermining the expansion as 
was used in iBietenholzl l|2006t ) and iBietenholz et al.1 (|2001bl . 
Il99ll ). and repeat a brief description here for the convenience 
of the reader. The goal is to determine the overall or average 
expansion rate of the radio nebula. Since there are few well- 
defined compact features, the expansion is measured not by 
determining the proper motion of individual features, but by 
determining an overall scaling between a pair of images by 
least-squares . This was accomplished by using the MIRIAD 
jSault et al.| [r995') task IMDIFFQ which determines how to 
make one image most closely resemble another, by calculat- 
ing unbiased estimators for the scaling in size, e, the scaling 
and the offset in flux density, A and h respectively, and the 
offsets in RA and deck, x and y respectively, by least squares. 
Our chief interest is in the expansion factor, e, but because 
of uncertainties in flux calibration, absolute position, and 
image zero-point offsets caused by missing short spacings, 

* As a check we also re-implemented the IMDIFF algorithm as 
an AIPS script. The results from this re-implementation, which 
uses different image-interpolation and minimisation and schemes, 
were consistent to within our uncertainties. 



all five parameters needed to be determined. This method 
was originally developed bv lTan fc GuTil l|l985h . 

In order not to contaminate our results on the expan- 
sion, we used an image made from the 2006 data set with- 
out the added D-configuration data from 1985. Although 
the structure at relatively small spatial scales, namely the 
filamentary structure visible in Fig. [TJ was well sampled by 
our A and B-configuration VLA observations, the largest 
scale structure was less well sampled, and may thus not 
have been sufficiently well recovered by the deconvolution. 
We will therefore obtain more reliable results on the ex- 
pansion if we use high-pass flltered images, which exclude 
the poorly sampled structure on the largest scales. To make 
the images for determining the expansion, we performed an 
initial high-pass filtering by using only visibility data at 
u-v distances > 6 KA. Then, we used CLEAN deconvolu- 
tion, since the positivity constraint applied in MEM decon- 
volution is not applicable to such high-pass filtered images. 
Finally, the resulting CLEAN images were further high-pass 
filtered in the image plane by applying a Gaussian high- 
pass filter of FWHM 15", and then the region exterior to 
G21. 5-0.9 was blanketfl. The 1991 1.5-GHz data were re- 
reduced in a manner consistent with the 2006 data, and a 
similarly high-pass filtered image was made. The FWHM 
resolution of this data set was 1.5" x 1.2" at —5°, the peak 
brightness was 6.2 mjy bm~^, and the rms background was 
0.12 mJy bm-^ Before running IMDIFF, the 2006 image 
was convolved to the resolution of the 1991 one. The fitting 
region encompassed ~2.8 square arcminutes or ~5000 beam 
areas. 

The fitted value of the expansion parameter, e, between 
epoch 1991.5 and 2006.3 was 1.017, in other words, the 
nebula has expanded by 1.7% over the 14.8 yr period be- 
tween the observations. The rms residual to the expansion 
fit was ~0.11 mJy bm~^, which is comparable to the com- 
bined background rms in the two images. 

The uncertainty of the expansion parameter is difficult 
to estimate reliably. The purely statistical uncertainty will 
be small due to the large number of image elements (~ 5000 
beam areas). However, the image elements are not strictly 
independent, so the expansion uncertainty is likely domi- 
nated by systematics such as deconvolution errors. We es- 
timate our uncertainty from the largest differences in the 
value of the expansion parameter, e, observed between IMD- 
IFF and the AIPS minimisation and between different lev- 
els of high-pass filtering. Our uncertainty is therefore not a 
statistical Icr value, but rather an estimated lo confidence 
interval. 

Since the two images are taken at different frequencies, 
caution is warranted when comparing them. However, as 
mentioned, we think it highly unlikely that any possible 
spectral index variations might cause a spurious apparent 
change in size with frequency for such a complex source 
structure. We therefore believe our measured expansion fac- 
tor accurately reflects true spatial expansion of the nebula. 
Expansion by a factor of 1.7 ± 0.3% over a 14.8 yr period 
suggests an expansion rate of 0.11 ± 0.02 % yr"'^, and in 
the absence of any acceleration or deceleration, an age of 



^ We use a slightly modified version of IMDIFF which treats 
blanked regions correctly 
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870_j5g yr. We note here that a small amount of accelera- 
tion, and therefore a slightly older age is expected. We dis- 
cuss this subject further in § [5] below. 

4.2 Spectral Index Image 

Having established the expansion rate, we can correct the 
1.5 GHz image for the expansion between 1991 and 2006, 
and then form a spectral index image with the two images, 
i.e., between 1.5 and 4.9 GHz. For the spectral index image, 
we used different images than for the expansion which are 
filtered neither in the visibility or the image-plane. The 5- 
GHz image was made incorporating the small amount of 
D-configuration data from 1985, after scaling of the latter, 
as described in § [2] above in order to recover the large-scale 
structure as accurately as possible. We show our spectral 
index image in Fig. O The average spectral index over our 
image (with spectral index, a, defined so that Sv oc i/") 
wa;0 +0.08to:S9- The rms of a over the image was 0.14. The 
average valu e of a is consistent w ith the earli er integrated 
value s of l|Morsi fc Reichlll987l ) and -0.02 (jSaher et all 
1 19891 ). 

The spectral index is relatively uniform over the neb- 
ula. Although a fairly large range of spectral index is visible 
in the image, a large uncertainty in the spectral index is 
unavoidable in regions of relatively low signal-to-noise. The 
formal uncertainty in the spectral index is, on average, 0.1 
over the central ~40" of the nebula, and as there are many 
image elements, a number of deviations exceeding 2g, i.e., 
0.2 in Of, are to be expected merely because of the noise. The 
regions showing larger deviations from the mean spectral in- 
dex are near the edge of the nebula where the uncertainty is 
higher. In addition, a wide band of slightly steeper spectral 
index (Aq ~ 0.15) running almost east- west is seen across 
the center of the nebula. This "feature" is almost certainly 
an artefact due to slight deconvolution errors due to incom- 
plete sampling near the center of the u-v plane. 

No systematic deviations from the mean spectral index 
are visible near the edge of the nebula, with a the mean 
spectral index of a strip of width ~10" around the circum- 
ference of the nebula differing only insignificantly (by —0.03) 
from the overall mean. No strong deviations are seen near 
the filamentary structures. The region near the pulsar has 
an a larger (i.e., the spectrum is more inverted) than that 
of its surroundings by ~ -1-0.1 and the brightest region has 
a slightly smaller a (i.e., the spectrum is slightly steeper). 
Although both of these deviations may be real, they are not 
larger than the uncertainties due to the noise and deconvo- 
lution errors. 



5 DISCUSSION 

We have made new, high-resolution 5-GHz radio images of 
the pulsar wind nebula, G21.5— 0.9, whose pulsar was only 
recently detected and found to have a very high spindown 
luminosity. By comparing our image to an earlier image at 

^ The uncertainty is the estimated systematic uncertainty due to 
flux calibration and a possibly slightly incomplete recovery of the 
total flux density in the 1.5 GHz image 



1.5 GHz we determined both the expansion rate and the 
radio spectral index of the nebula . 

Our image clearly shows that G21.5— 0.9 has a filamen- 
tary structure in the radio, similar to that seen in radio 
images of the Crab and of 3C 58. Large loop-like struc- 
tures are seen especially towards the northeast and south- 
west. Such s t ructu re seems to be present in all young PWNe. 
ISlane et al.l (|2004l ) suggested that they were loops torn by 
the kink instability from the toroidal field produced by the 
pulsar. 

The nebula is expanding fairly rapidly: our estimate 
of the expansion rate between 1991 and 2006 is 0.11 ± 
0.02% yr~^. At 5 kpc this rate corresponds to a speed (wrt. 
the centre) of 910 ± 160 km s~^. This is similar to although 
a little slower than the expansion speeds seen for two other 
young PWNe: the Crab Nebula and 3C 58. 

If we assume constant velocity expansion, then the age 
of G21.5-0.9 is 870ti^[; yr. This is considerably younger 
than the characteristic age of the pulsar, which is 4800 yr. 
In fact, the nominal age of 870 yr is ~100 yr younger then the 
Crab pulsar and would make PSR J1833— 1034 the youngest 
known Galactic pulsar with the possible exception of PSR 
J1846— 0258 in the remnant Kes 7fJl|. We note, however, that 
the expansion of a young PWN is expected to accelerate 
slightly as it expands into the freely expanding ejecta for 
the first -^1000 yr of its hfe (|van der Swaluw et al.1 1200 ll : 
iBlondin et a"l]|200ll : [junlll998l ). Such acceleration is in fact 
observed in the case of the Crab Nebula (|Trimbld [l968l : 
lNugent|[l99i : iBietenholT et al.|[T99ll ). If G21.5— 0.9 is simi- 
larly accelerated, then its age would be somewhat older than 
870 yr. In fact, it would have very nearly the same age as 
the Crab nebula, since both the constant-v elocity expansion 
age o f both the Crab's synchrotron nebula (|Bietenholz et all 
I1991I ) is almost the same as that we have determined for 
G21.5— 0.9. The only other pulsar known to be younger 
than this is that in 0540 — 69, whose age is also deter- 
mined by measuring the expansion speed of the nebula, al- 
beit from the optical spec trum rather than from radio imag- 
ing (iKirshner et al.ll 19891 ). With the exception of the pulsar 
spindown age, which often differs from the true age of the 
pulsar by factors of a few, this is the first age measurement 
for G21. 5-0.9. 

The supernova responsible for creating PSR 
J1833-1034 and G21. 5-0.9 would then have exploded 
around the year 1140 AD. It might be considered surprising 
that no historical supernova was reco rded at tha t time 
near G2 1.5 -0.9, ahhough we note that IWang et all (119861 ) 
suggested that G21.5— 0.9 might be the remnant of the 



No expansion-speed measurement is available for Kes 75. 
The spindown age o f PSR J1846-0258 is 723 yr (e.g., 
iLivingstone et aU I2OO6I ) . Arguments have been advanced that 
Kes 75, and thus its pulsar, are considerably older. These argu- 
ments, howev er, were based in part on a large distance to Kes 75. 
Very recently, iLeahv &: TianI ||2007|) , have shown that a consider- 
ably smaller distance is likely for PSR J1846— 0258, and thus an 
age nearer the spindown age is quite possible. The other possi- 
ble exception is PSR J0205-I-6449 in 3C 58, which has often been 
associated with the supernova of 1181 A.D., which would give 
it an age of ~825 yr. A number of arguments, however, suggest 
that 3C 58 and PSR J0205+6449 are conside rably older (se c for 
example the discussions in Bietenholz 2006 or IChevalierll2005l) 
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Figure 3. An image of the radio spectral index of G21.5— 0.9 between 1.5 and 5 GHz. The spectral index (Si, oc u") is plotted in 
greyscale, and regions where the uncertainty in the spectral index were larger than 0.2 are blanked. The FWHM resolution is 1.5 X 1.2" 
at p. a. —5°. The two contours show the 5 and 10% contours of the image at 5 GHz for reference (note that this image is at lower 
resolution than the one in 

Fig-CJ- The cross again denotes the pulsar position. 



"guest star" (supernova) of 48 BC, making it ~2100 yr old, 
which is compatible with our measurement only if strong 
acceleration has occurred since the supernov a explosion. 
However, as pointed out bv lCamilo et al.l l|2006l '). the hydro- 
gen column density i mplied by various X-r a y measurements 
(~ 2 X 10^^ cm-^ [Warwick etall 1200 ll : ISafi-Harb etall 
|2001^ in the direction of G21.5— 0.9 corresponds to sufflcient 
visual extinction (10-11 mag) for it to be unlikely that the 
supernova event was noticed. 

We also produced a radio spectral index image, which 
shows that G21.5— 0.9's spectrum is fairly uniform, with 
an average spectral index between 1.5 and 5 GHz of a = 
-|-0.08to:o9- G21. 5-0.9 is similar in this regard to two other 
young PWNe, the Crab and 3C 58, whose radio spectra 
have been found to be remarkably uniform over the nebu- 
lae. The spectral uniformity suggests there is a single source 
for the relativistic electrons responsible for the radio emis- 
sion, namely the pulsar. 



One might also expect some radio emission from elec- 
trons accelerated in the shocks produced by the interaction 
of the expanding supernova ejecta with their surroundings, 
and those bounding the pulsar nebula as it expands into the 
ejecta. Such shocks are expected to give rise to radio spec- 
tra steeper than those observed in PWNe (typically with 
Q < —0.5). No such regions are seen in G21.5— 0.9. In par- 
ticular, on our image, a region of width ~10" around the 
edge of the nebula has a mean spectral index which differed 
by < 0.03 from the mean over the entire nebula. 

As mentioned earlier, a faint "halo" of ~2.5' in di- 
ameter is seen surrounding the PWN in the X-ray, with 
abou t ~1% of the la ttcr's surface brightness ("Bocchino ct ajj 
20051: iMatheson fc S afi-Harb 200_5; S afi-H arb ct al. 200l|; 



Slane et al.l 



2000l ) . This halo might be X-ray emission 
associated with the o uter shock, although it s nat ure is 
not yet clear. Both IMatheson" fc Safi-Harbl l|2005l ) and 
iBocchino et al.l l|2005l ) argue that most of the halo X-ray 
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emission is not from the outer shock, but is due to dust 
scattering. 

Is there any radio emission associated with this halo? 
Previous radio observations have not shown the presence of 
a corresponding radio halo to a la surface brightness limit 
of 4 X 10"^^ W m~^ H z ~^ sr~^ (.^30 mj y arcmin-^) at 
1 GHz (jBock et alJboOll : ISlane et al.ll2000h . We produced 
a wide-field radio image from our new radio data to search 
for such extended emission, but found no features above the 
noise outside the region shown in Fig. [1] We note, how- 
ever, that our image background rms of only about 1% of 
the brightness of the PWN is probably not low enough to 
detect a halo comparable to the X-ray one, and is in fact 
somewhat higher than existing limits on the halo radio emis- 
sion. Deeper, low-frequency radio observations of G21.5— 0.9 
should be undertaken to determine, or set lower limits to, 

the radio brightness of the halo. 

It has already been poin ted out, by IWoltier et al.l 

|l993) and ISalvati et al.l (|l998ll . that interpreting the rel- 
atively low break frequency (~100 GHz) measured in the 
broadband spectrum of G21.5— 0.9 (see Boc k et al.l l200ll : 
iGallant fc TuSll998l : ISalter et al.l [l989,) as being caused 
by synchrotron ageing implies either an unrealistically high 
magnetic field or an unrealistically large age. For example, 
if we take the minimum ener gy field calculated from syn- 
chrotron theory of -^440 fiG f Slane et al . 2000), and take 
the break freque ncy to be 100 GH z, one would arrive at an 
age of 13,000 yr (iBock et al.ll200il ). Our measurement of the 
expansion speed now provides firm observational evidence 
against such large ages. A much larger magnetic field, on 
the other hand, is not feasible on energetic grounds. Pos- 
sibly t he break in the spectrum is not due to synchrotron 
ageing. [Fleishman fc Bietenholj |2003) showed that, in the 
case of a tangled magnetic field, several breaks occur in the 
synchrotron spectrum even for a single power-law distribu- 
tion of relativistic electrons, and these breaks can naturally 
explain the broadband spectra of PWNe, both ones like the 
Crab with break frequencies in the infrared, and ones like 
G21.5— 0.9 or 3C 58, which have much lower break frequen- 
cies. Such a model might be applicable to G21.5— 0.9 as well. 

Can we identify any radio features associated with the 
termination shock in the pulsar wind? The termination 
shock's radius can be estimated by balancing the ram pres- 
sure of the pulsar wind against the nebular pressure, and 
one finds that the radius of the termination shock is ^l'.'6 
rj where t? is the fra ctional solid angle covered by the wind 
ifCamilo et al.ll2006l . see also Slane et al. 2000, who obtained 
a termination shock angular radius of l'.'5 rj for input pa- 
rameters somewhat difi'erent than ours.) . 

The pulsar in G21.5— 0.9 is surrounded by a region of 
somewhat lower radio surface brightness, and the possible 
elliptical structure mentioned in § [3] (see Figure [2]) , which 
has a radius of ~4". In high-resolution X-ray images, there 
is a central, bright region which has an elliptical shape, elon- 
gated al ong the northwest-s outheast direction, with a size of 
7" X 5" (jCamilo et al.l|2006l , see also Matheson & Safi-Harb 
2005; Slane et al. 2000). If both these radio and X-ray struc- 
tures are interpreted as being associated with the pulsar 
wind termination shock, then a radius for the latter of ~4" 
is suggested, which implies that 77 ~ 0.4. Given that the only 
two well-resolved termination shock structures, being those 
in the Crab Nebula and in Vela, seem to lie predominately 



in the pulsar's equatorial plane, and thus have rj < 1, it 
seems this interpretation is not unreasonable. A more con- 
clusive interpretation will have to await higher-resolution 
radio and/or X-ray observations. 
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